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Abstract 

Aims. We present an extension to our XMM-Newton X-ray source catalogue of M 31 containing 39 newly found sources. In order to 
classify and identify more of the sources we search for X-ray time variability in XMM-Newton archival data of the M 31 centre field. 
Methods. As a source list we used our extended catalogue based on observations covering the time span from June 2000 to July 
2004. We then determined the flux or at least an upper limit at the source positions for each observation. Deriving the flux ratios for 
the different observations and searching for the maximum flux difference we determined variability factors. We also calculated the 
significance of the flux ratios. 

Results. Using hardness ratios, X-ray variability and cross correlations with catalogues in the X-ray, optical, infrared and radio 
regimes, we detected three super soft source candidates, one supernova remnant and six supernova remnant candidates, one globular 
cluster candidate, three X-ray binaries and four X-ray binary candidates. Additionally we identified one foreground star candidate 
and classified fifteen sources with hard spectra, which may either be X-ray binaries or Crab-like supernova remnants in M 31 or 
background active galactic nuclei. The remaining five sources stay unidentified or without classification. Based on the time variability 
results we suggest six sources, which were formerly classified as "hard", to be X-ray binary candidates. The classification of one other 
source (XMMM31 J004236.7-I-41 1349) as a supernova remnant, has to be rejected due to the distinct time variability we found. We 
now classify this source as a foreground star. 
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1. Introduction 



An ideal target for a search for time variability of X-ray sources 
is t he bright Local Group spiral galaxy M 3 1 (distance 780 
kpc. lHollandll998l:IStanek & Garnavichll998h with its moderate 
Galac tic foreground absorption (A^h= 7x10^*^ cm"^, I Stark et all 
Il992h . 

The Einstein X-ray observatory found 16 sources in M 31, 
which showed v ariability comparing the individual observations 
with each other (van Spevbr oeck et al.l 19791: ICollura et al.ll990t 
iTrinchieri & Fabbiano 1 99 ll hereafter TF91). iMmimetalJ 
( I1993L hereafter PFJ93) compared ROSAT HRI to previous 
' Einstein observations and found several variable sources. The 
two ROSAT PS PC surveys of M 3 1 , covered the entire galaxy 
and were separated by about one vear. lSupper et al.l(ll997il20"0]1) 
found, that the intensity of 34 sources varied significantly be- 
tween the observations. 

iGarcia et al.l ( l2000l) reported on first observations of the 
nuclear region of M 31 with Chandra. They found that the 
nuclear source has an unusual X-ray spectrum compared to 
the other point sources in M 31. Source catalogues, based on 
Chandra observation s , of th e ce ntral fie l d of M 31 are pro- 
vided by iKong et al.l (|2002|) and iKaarej (|2002|) . Three differ- 
ent M 31 disk fields, spanning a range of stellar populations, 
were observed by Chandra to compare their point source lu- 
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minosi ty functions to that of the galaxy's bulge dKong et al.l 
12003 al). In a synoptic study of the disk 0.9 square degree) of 
M 3 1 . IWilliams et al.l (l2004l) measured the mean flux and long- 
term light curves for 166 objects. At least 25% of the sources 
show significant variability. Bright X-ray binaries (XRBs) in 
globular clusters and supersoft sourc es (SSSs) and quasisoft 
sources (Q SSs) were investigat ed by iDi Stefano et a U (120021 
|2004 and iGreiner et all (|2004. Th e disco v ery of an X-ray 
nova w as reported bv IWilliams et all ( I2005ah . IVoss & Gilfanovl 
(120071) used Chandra data to examine the low mass X-ray bi- 
naries (LMXBs) in the bulge of M 31. Good candidates for 
LMXBs are the so-called transient sources. Studies of tran- 
sient sources_Jn_M 31 are presen ted in nu merous papers, 
e. g. IWilliams etall (l200 6b), Trudol vubov et al. (2006a here- 
after TPC06), IWiUiams et al.l ( l2005bl) . IWilUams et al.1 (l2006al 
hereafter WGM06). Using XMM- Newton and Chandra data, 
iTrudolvubov & Priedhorskvl (l2004h detected 43 X-ray sources 
coincident with globular cluster candidates from various opti- 
cal surveys. They studied their sp ectral properti e s, time vari- 
ability and luminosity functions. lOsborne et al.l ( 1200 lb used 
XMM-Newton Performance Verification observations to study 
the variability of X-ray sources in the central region of M 31. 
They found 116 sources brighter than a limiting luminosity of 
6x10-'^ ergs"' and examined the ~ 60 brightest sources for 
periodic and non-periodic variability. At least 15% of these 
sources appea r to be v ariable on a time scale of several months. 
iBarnard et aTl (l2003al) used XMM-Newton to study the X-ray 
binary RX J0042.6+4115 and suggested it as a Z-source. lOriol 
(,2006i) s tudied the population of S SSs and QSSs with XMM- 
Newton. |^udolyubov|eral](|20063) provide a study of the bright 
sources in the central region of M 31, including spectral proper- 
ties, variability and source classification. It is based on the same 
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XMM-Newton observations analy sed in this paper. Recently 
iTrudolvubov & Pri edhorskvl (|2007|) reported the discovery of 
21 7s pulsations i n a bright persistent SSS. 

iPietsch et afl (l2005bi hereafter PFH2005) prepared a cata- 
logue of M 3 1 point-like X-ray sources analysing all observa- 
tions available at that time in the XMM-Newton archive which 
overlap at least in part with the optical D25 extent of the galaxy. 
In total, they detected 856 sources. The centre part of the galaxy 
was covered four times with a separation of the observations 
of about half a year starting in June 2000 (some of the re- 
gions at the boundary of the centre area are even covered five 
times). PFH2005 give only source properties derived from an 
an alysis of the comb i ned ce ntre observations. In follow-up work 
m IPietsch &Haberll(l2005h searched for X-ray burst sources in 
globular cluster (GIC) sources and candidates and i dentified two 
X-ray b ursters and a few more candidates, and (ii) lPietsch et alj 
(l2005ai hereafter PFF2005) searched for correlations with op- 
tical novae. They identified seven SSSs and one symbiotic star 
from the list of PFH2005 with optical novae and identified one 
additional XMM-N ewton s ource with an optical nova. This work 
was continued in Pi etsch et al. (200 7, hereafter PHS2 007). 

Similar to the M 33 work of IPietsch et al.l (|2004|) PFH2005 
used the hardness ratios, i. e. X-ray colours, and correlations 
with sources in other wav elength regimes t o iden tify and clas- 
sify the detected sources. iMisanovic et al.l ( l2006l) showed, for 
a source population study of M 33, that X-ray flux variability 
on different time scales allows us to further distinguish between 
different source classes. Phenomena such as bursts of X-ray bi- 
naries, flares of stars or the periodic variability of pulsars oc- 
cur on short time scales and can therefore be observed during 
one single observation. On the other hand there is long term 
variability. Those time scales can be covered, comparing dif- 
ferent observations of the same source. In the field of view of 
M 3 1 there are mainly two source classes, which are known to 
show strong variability (variability factor > 10) on time scales 
of years. These are X-ray binaries and SSS. Among the active 
galactic nuclei (AGNs) narrow-line Seyfert 1 galaxies and BL 
Lac objects show the strongest variabilities. However only a few 
narrow-line Seyfert 1 galaxies are known in the entire sky, which 
show flux variability factors larger than 10 on time scales of half 
a year up to several years. Hence it is very unlikely that one of 
the strongly variable sources in M 3 1 would be an AGN. 

In this paper we report a search for new X-ray sources in 
the XMM-Newton observations to the M 31 centre to extend 
the source catalogue of PFH2005 and a time variability analysis 
of all the M 31 centre sources. In Sect. 2 information about the 
used observations and accomplished analysis is provided. Sect. 3 
describes the source catalogue extension. The results of the tem- 
poral variability analysis are discussed in Sect. 4. Discussion of 
the individual source classes, including X-ray identifications, are 
provided in Sect. 5. We draw our conclusions in Sect. 6. 

2. Observations and analysis 

For our analysis we used the archival XMM-Newton observa- 
tions of the central region of M 31, obtained from June 2000 
to January 2002 (from observations si and nl only sources 
which lie in the intersection with at least one of the other 
observations are included). In addition, we analysed the July 
2004 monitoring observations of the low mass X-ray binary 
RX J0042.6+4115 (PI Barnard), which are pointed 1.1' to the 
west of the M 31 nucleus position. Thus we have a time span 
of about four years for examining time variability. Details of 
the observations can be found in Table[T]which shows the M 31 



field name (Col. 1), the identification number (2) and date (3) 
of the observation and the pointing direction (4, 5). Column 6 
contains the systematic position offset. For each EPIC camera 
the used filter and the exposure time after screening for high 
background is given. To achieve comparable images and results 
we adapted the same background screening as in PFH2005 
for the newly added observations. We had to reject ObsID 
0202230301, because it shows high background throughout the 
observation. To increase the detection sensitivity we merged 
the data of ObsID 0202230201, 0202230401 and 0202230501 
after correction of the position offsets. The combination of these 
three observations is called "b". 

We searched for sources in "b", which were not visible in 
the X-ray wavelength regime about 2.5 years earlier. In addition 
we reexamine observations cl, c2, c3 and c4 individually, to 
search for sources not included in the PFH2005 catalogue, 
which - besides source 856 - was based on an analysis of the 
merged images of observations cl to c4. 

The data analysis was performed using tools in the XMM- 
Newton Science Analysis System (SAS) v6.6.0 and v7.0.0, 
EXSAS/MIDAS 03OCTJEXP and FTOOLS v6.0.6 software 
packages, the imaging application DS9 v4.0b7 together with the 
funtools package, the mission count rate simulator WebPIMMS 
v3.6a and the spectral analysis software XSPEC vll.3.1. 

2.1. Images 

We used five energy bands: (0.2-0.5) keV, (0.5-1.0) keV, (1.0- 
2.0) keV, (2.0^.5) keV, and (4.5-12) keV, to create images, 
background images and exposure maps for PN, MOSl and 
MOS2 and masked them for acceptable detector area. For PN, 
the background maps contain the contribution from the "out 
of time (OOT)" events (parameter withootset=true in task 
esplinemap). 

Figure [T] shows logarithmically scaled XMM-Newton EPIC 
low background images of the M 31 centre observations inte- 
grated in l"xl" pixels combining data from the PN, MOSl and 
MOS2 cameras in the (0.2-4.5) keV XID band. The data are 
smoothed with a 2D-Gaussian of FWHM 5", which corresponds 
to the point spread function in the centre of the field of view 
(FOV). Figure |2] gives a zoom-in to the crowded centre region. 

2.2. Source detection 

We searched for sources using simultaneously 5x3 images 
(5 energy bands for each EPIC camera). A preliminary source 
list created with the task eboxdetect with a low likelihood 
threshold (likemin = 5) was used as starting point for the 
task emldetect (v. 4.44.19). We used parameters nmaxfit = 
2 and fitextent = true. The parameter extentmodel was 
set to beta and we only allowed multi-PSF fitting if the likeli- 
hood was larger than 10. Setting the parameter withtwostage 
to true the program checked in a second run whether split- 
ting extended sources into two point-like sources (nmulsou = 
2) would achieve a more reliable result. 

For most sources, band 5 just adds noise to the total count 
rate. If converted to fluxes this noise often dominates the total 
flux due to the high ECF. To avoid this problem we calculated 
count rates and fluxes for detected sources in the "XID" (0.2- 
4.5) keV band (bands 1 to 4 combined). While for most sources 
this is a good solution, for extremely hard or soft sources there 
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Table 1. XMM-Newton log of archival M 31 observation overlapping with the optical Dos ellipse (proposal numbers 010927, 
011257 and 015158). 



M 31 field 


Obs. id. 


Obs. dates 


Pointing 


direction 


Ofl^set • 


EPIC PN 


EPIC MOSl 


EPIC M0S2 








RA/Dec 






Filter^ 


j-f 

i e xp 


Filter^ 


i exp 


Filter^ 


^ exp 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Centre 1 (cl) 


0112570401 


2000-06-25 


0:42:36.2 


41:16:58 


-1.9, +0.1 


medium 


26.40 


medium 


29.92 


medium 


29.91 


Centre 2 (c2) 


0112570601 


2000-12-28 


0:42:49.8 


41:14:37 


-2.1, +0.2 


medium 


9.81 


medium 


12.24 


medium 


12.24 


Centre 3 (c3) 


0112570701 


2001-06-29 


0:42:36.3 


41:16:54 


-3.2,-1.7 


medium 


27.65 


medium 


27.65 


medium 


27.65 


North 1 (nl) 


0109270701 


2002-01-05 


0:44:08.2 


41:34:56 


-0.3, +0.7 


medium 


54.78 


medium 


57.31 


medium 


57.30 


Centre 4 (c4) 


0112570101 


2002-01-06/07 


0:42:50.4 


41:14:46 


-1.0, -0.8 


thin 


60.79 


thin 


63.31 


thin 


63.32 


South 1 (si) 


0112570201 


2002-01-12/13 


0:41:32.7 


40:54:38 


-2.1,-1.7 


thin 


53.45 


thin 


53.76 


thin 


53.73 


(bl)* 


0202230201 


2004-07-16 


0:42:38.6 


41:16:04 


-1.3,-1.2 


medium 


18.30 


medium 


19.40 


medium 


19.40 


(b2) 


0202230301 


2004-07-17 


0:42:38.6 


41:16:04 


-1.0, -0.9 


medium 


0.0 


medium 


0.0 


medium 


0.0 


(b3)* 


0202230401 


2004-07-18 


0:42:38.6 


41:16:04 


-1.7,-1.5 


medium 


13.80 


medium 


17.90 


medium 


17.90 


(b4)t 


0202230501 


2004-07-19 


0:42:38.6 


41:16:04 


-1.4,-1.8 


medium 


8.90 


medium 


10.20 


medium 


10.20 



Notes: 

* : Systematic offset in RA and Dec in arcsec determined from conelations with 2MASS, USNO-B 1 and Chandra catalogues 
^ : All observations in full frame imaging mode 

^ : Exposure time in units of ks after screening for high background used for detection 
^ : Combination of the three observations is called b (see text) 



may still be bands just adding noise. This then may lead to rate 
and flux errors that seem to wrongly indicate a lower source 
significance. A similar effect occurs for the all instrument rates 
and fluxes if a source is mainly detected in one instrument (e.g. 
soft sources in PN). 

We accepted sources which have a likelihood above 6 in the 
combined fit. We rejected spurious detections in the vicinity of 
bright sources. In regions with a highly structured background 
the SAS detection task emldetect registered some extended 
sources. We also rejected these "sources" as spurious detections. 
Two sources were added manually: Source 87 1 was first detected 
as nova WeCA PP-N2001-12 and in flie POINT - AG APE variable 
star catalogue ( I An et al.ll2004.lAn et alj (|2004|) pro pose the hard 
X-ray transient [OBT2001] 3 (lOsborne et al.ll200ll) as a counter- 
part, which is source 287 in the PFH2005 catalogue. PFF2005 
showed that several points speak against this identification and 
that a faint SSS close to the position of [PFH2005] 287, which 
is only visible during observation c4, is a more reliable counter- 
part. Source 885 (see Table 3) is clearly visible in observation b 
(see Fig. O and we could not find any reason, why emldetect 
did not automa tically find it. As t he source was already reported 
with Chandra jKong et al.ll2002L r2-41), we took the Chandra 
position to derive the source parameters, using emldetect with 
fixed position. 

Our source catalogue extension only contains sources not al- 
ready found by PFH2005. These sources were ordered according 
to increasing right ascension for each observation individually. 
Finally we merged the source lists and numbered the sources 
consecutively. If a source was detected in more than one obser- 
vation, we took the source parameters from the first observation, 
in which it was detected. As this catalogue is an extension of the 
source Ust of PFH2005, new sources start with number 857. 

To classify the source spectra we computed four hardness ra- 
tios from the source count rates. These hardness ratios and errors 
are defined as 

HRi = — '- and EHRi = 2-^-^ — —, (1) 

B/+1 + Bi (Bm + B,)2 

for ; = 1 to 4, where B, and EBj denotes count rates and cor- 
responding errors in band ; as defined above. The identification 



and classification criteria are given in Table |2] The source cata- 
logue extension is presented in Sect. |3](see Table 3). 



2.3. Variability calculation 

To examine the time variability of each source, we determined 
the XID flux - or at least an upper limit for the XID flux - 
at the source position in each observation . We used the task 
emldetect (v. 4.60) setting the parameter xidfixed = yes 
which forced emldetect not to alter the source positions in cal- 
culating the total flux. To obtain fluxes and upper limits for all 
sources in our input list we set the detection likelihood threshold 
to 0. 

Merging the source catalogue extension (see Sect. 3) with 
the source catalogue of PFH2005 we generated a starting list 
for our variability analyses. This starting list only contains the 
number and position of each source. To give correct results the 
task emldetect has to process the sources from the brightest to 
the faintest one. Therefore we first had to order the sources in 
each observation by detection likelihood. For sources not visible 
in the observation we arbitrarily set the detection likelihood to 1 . 
This list was used as input for a first emldetect run. This way 
we achieved an output list, in which a detection likelihood was 
allocated to every source. To finally examine the sources ordered 
by detection likelihood, a second emldetect run was necessary. 

We only accepted XID fluxes, which are at least three times 
larger than their 1 cr errors. Otherwise the triplicated error was 
used as an upper limit. The largest XID flux of each source was 
derived, excluding upper limit values (column fmax in Table 5). 
Comparing the XID fluxes of the diff'erent observations with 
each other, we calculated the significance of the difference 



+ cr'- 



(2) 



(column svar in Table 5) and the ratio of the XID fluxes Fvar = 
Fmax/Fmin (column fvar in Table 5), if Fmax was not an upper 
limit. Fmax and Fmin are the source XID maximum and mini- 
mum (or upper limit) flux and iXmax and cr^i„ are the errors of 
the maximum and minimum flux, respectively. The results of the 
time variability analyses are discussed in Sect. 4. 
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Figure 1. Logarithmically scaled XMM-Newton EPIC low background images integrated in 1" pixels of the M 31 centre observa- 
tions combining PN and MOSl and MOS2 cameras in the (0.2-4.5) keV XID band. The data are smoothed with a 2D-Gaussian of 
FWHM 5", which corresponds to the point spread function in the centre area. The images are corrected for unvignetted exposures. 
Contours are at (2,4, 8, 16, 32) x 10"* ct s"' pix"' including a factor of two smoothing. Sources from the combined catalogue are 
marked in the outer area. The inner area is shown in detail in Fig.|2l 



3. Source catalogue 

PFH2005 reported 265 sources in the centre field of M 31. Our 
catalogue extension contains 39 sources. Four are detected in 
observation cl, eight in observation c3, thirteen in c4 and twenty 
one in "b". 

The source parameters are summarised in Table 3 (EPIC 
combined products and products for EPIC PN, MOSl and 
MOS2, separately). 



With the exception of the newly added XMM-Newton source 
name (column 77, see below) Table 3 is structured in the same 
way as Table 2 from PFH2005. It gives the source number 
(Col. 1), detection field, from which the source was entered into 
the catalogue extension (2), source position (3 to 9) with Icr un- 
certainty radius (10), likelihood of existence (11), integrated PN, 
MOS 1 and MOS2 count rate and error (12,13) and flux and error 
(14,15) in the (0.2^.5) keV XID band, and hardness ratios and 
errors (16-23). Hardness ratios are calculated only for sources 
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Figure 1. (continued) Logarithmically scaled XMM-Newton EPIC low background images integrated in 1 "pixels of the M 3 1 centre 
observations combining PN and MOSl and MOS2 cameras in the (0.2-4.5) keV XID band. 



for which at least one of the two band count rates has a signif- 
icance greater than 2cr. Errors are the properly combined sta- 
tistical errors in each band and can extend beyond the range of 
allowed values of hardness ratios as defined previously (-1.0 to 
1.0). The EPIC instruments contributing to the source detection, 
are indicated in the "Val" parameter (Col. 24, first character for 
PN, second MOSl, third MOS2) as "T", if inside the FOV, or 
"F", if outside of FOV. There are 8 sources at the periphery of the 
FOV where only part of the EPIC instruments contribute. The 
positional error (10) does not include intrinsic systematic errors 
which amount to (J! 5 (see PFH2005) and should be quadratically 
added to the statistical errors. 

Table 3 then gives for EPIC PN, exposure (25), source ex- 
istence likelihood (26), count rate and error (27, 28) and flux 
and error (29, 30) in the (0.2-4.5) keV XID band, and hardness 
ratios and errors (31-38). Columns 39 to 52 and 53 to 66 give 
the same information corresponding to Cols. 25 to 38, but for 
the EPIC MOSl and MOS2 instruments. Hardness ratios for the 
individual instruments were again screened as described above. 
From the comparison of the hardness ratios derived from inte- 
grated PN, MOSl and MOS2 count rates (Cols. 16-23) and the 
hardness ratios of the individual instruments (Cols. 31-38, 45- 
52 and 59-66) it is clear that combining the instrument count rate 
information yielded significantly more hardness ratios above the 
chosen significance threshold. 

Column 67 shows cross correlations with M 3 1 X-ray cata- 
logues in the literature. 

Our catalogue extension contains 23 until now unknown X- 
ray sources in M 31. The discussion of the results of the cross 
correlation is in Sect.|5] 

In the remaining columns of Table 3, we give cross correla- 
tion information with sources in other wavelength ranges. 



To identify the X-ray sources in the M 3 1 field we searched 
for correlations around the X-ray source positions within the 3 cr 
total X-ray error in the SIMBAD and NED archives and within 
several catalogues. In columns 68 to 73 of Table 3, we give ex- 
traction information from the USNO-Bl catalogue (name, num- 
ber of objects within search area, distance, B2, R2 and I magni- 
tude of the brightest object). To improve the reliability of identi- 
fications we used the B and R magnitudes to calculate 

log(:^) = log(/x) + (mB2 + mu2)/{2 X 2.5) + 5.37, (3) 

lopt 

following lMaccacaro et al.l (Il988i see Col. 74). 

The X-ray sources in the catalogue extension are identi- 
fied or classified based on properties in the X-ray (HRs, vari- 
ability) and of correlated objects in other wavelength regimes 
(Table 3, Cols. 75, 76). The criteria are summarised in Table |2] 
Identification and classification criteria are discussed in detail in 
Sect. 6 of PFH2005. As we have no clear hardness ratio criteria 
to select XRBs, Crab-like supernova remnants (SNRs) or AGNs 
we introduced a class <hard> for those sources. If such a source 
shows strong variability (i.e. V > 10) on the examined time 
scales it is likely to be an XRB. Fifteen sources are classified as 
<hard>. Five sources remain unidentified or without classifica- 
tion. 

The last column (77) of Table 3 contains the XMM-Newton 
source name as registered to the lAU Registry. Source names 
consist of the acronym XMMM31 and the source position as 
follows: XMMM31 Jhhmmss.s+ddmmss, where the right ascen- 
sion is given in hours (hh), minutes (mm) and seconds (ss.s) 
truncated to decimal seconds and the declination is given in de- 
grees (dd), arc minutes (mm) and arc seconds (ss) truncated to 
arc seconds, for equinox 2000. 

Only two sources from our catalogue extension (869, 863) 
are found as extended sources (see Table|4]and Sect.|5]). 
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Figure 2. Inner area of M 31 enlarged from Fig.[T] Contours are at (4, 8,16, 32, 64, 128, 256) x 10"* ct s ' pix ' including a factor of 
two smoothing. Sources from the combined catalogue are marked as 30"x30" squares. The images are ordered as follows: Centre 
1 (upper left). Centre 2 (upper right). Centre 3 (middle left). Centre 4 (middle right) and Centre B (lower left). 
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Source type"' Selection criteria Identified Classified 

fgStar log{4^) < -l.O md HR2-EHR2< 0.3 and HR3-EHR3< -OA or not defined 1 

/opl 

AGN Radio source and not classification as SNR from HR2 or optical/radio 

Gal optical id with galaxy 

GCl X-ray extent and/or spectrum 

SSS HRl < 0.0, HR2 - EHR2 < -0.99 or HR2 not defined, HR3, HR4 not defined 3 

SNR HRl > -0. 1 and HR2 < -0.2 and not a fg Star, or id with optical/radio SNR 1 6 

GIC optical id 1 

XRB optical id or X-ray variability 3 4 

hard HR2 - EHR2 > -0.2 or only HR3 and/or HR4 defined, and no other classification 15 



Notes: 

t : fg Star: foreground star, AGN: active galactic nucleus, Gal: galaxy, GCl: galaxy cluster, SSS: supersoft source, SNR: supernova remnant, GIC: 
globular cluster, XRB: X-ray binary 



Table 4. Extension properties of sources 863 and 869 





Source 


Extent 
arcsec* 


Ext. err."' 
arcsec* 


MELH* 


863 
869 


6.71 
6.39 


2.14 
1.12 


4.70 
5.05 



Notes: 

' : Extent error 

* : Maximum extent likelihood 

* : 1" corresponds to 3.8 pc at the assumed distance of M 31 

4. Variability 

Table 5 contains all information necessary to examine time vari- 
ability. The sources are taken from the combined catalogue 
(i. e. PFH2005 and Sect. |3]l. Sources are only included in the ta- 
ble, if they are in the FOV for at least two observations. Column 
1 gives the source number Columns 2 and 3 contain the flux 
and error in the (0.2-4.5) keV XID band. The hardness ratios 
and errors are given in columns 4 to 11. Column 12 shows cross 
correlations with M 31 X-ray catalogues in the literature. The 
next two columns contain the type of the source (13) and cross 
correlation information with sources in other wavelength ranges 
(14). The EPIC instruments contributing to the source detection 
in the cl observation, are indicated in the "cl_val" parameter 
(Col. 15, first character for FN, second MOSl, third MOS2) as 
"T", if inside the FOV, or "F", if outside FOV. Then the count 
rate and error (16,17) and flux and error (18,19) in the (0.2^.5) 
keV XID band, and hardness ratios and error (20-27) of the cl 
observation are given. Corresponding information is given for 
observation c2 (cols. 28-40), c3 (41-53), nl (54-66), c4 (67- 
79), si (80-92) and b (93-105). 

Column 106 indicates the number of observations in which 
the source is covered in the combined EPIC FOV. The maxima 
of the significance of variation and flux ratio (fvar_max) are 
given in columns 107 and 108. As described in Sect. 2.3 we 
only used detections with a significance greater 3 cr. Otherwise 
the 3 cr upper limit was used. Column 109 indicates the number 
of observations where we could only gain an upper limit. The 
maximum flux (fmax) and its error are given in columns 110 
and 1 1 1 . In a few cases we could not derive the maximum flux, 
because every observation only gives an upper limit. This can 
have two reasons: The first reason is that PFH2005 merged 



observations cl to c4 for source detection. Hence a faint source 
may not be detectable at the 3cr limit in the individual observa- 
tions. The second reason is, that in cases where the significance 
of detection was not much above the 3 cr limit, it can become 
smaller than the 3 cr limit when the source position is fixed. 
The source name, according to the lAU naming convention (see 
Sect. [3]), can be found in column 1 12. 

In Fig. [3] we plotted the variability factor (col. fvarjnax) 
of each source as function of its maximum flux (col. fmax) 
in the XID band. Identified sources are marked with big sym- 
bols, whereas classified sources are indicated by small symbols. 
Source numbers from PFH2005 and Sect. [3] are indicated for 
sources with flux variability above 5 or maximum flux above 
8x10"'^ erg cm"^ s"'. In this region only, variability can help 
distinguish between foreground stars or SNRs, or to decide if a 
source classified as hard is an AGN or a XRB. Sources with a 
statistical significance of the variability below 3 are marked in 
green (grey). 

Figure [3] clearly shows that most of the variable sources are 
XRBs or XRBs in GIC or candidates of these source types. In ad- 
dition there are a few SSS candidates, and even some SNR candi- 
dates showing pronounced temporal variability. These SNR can- 
didates are discussed in Sect. |5] as they should not show time 
variability. The sources classified or identified as AGNs or fore- 
ground stars all show Fvar < 4, besides the new foreground star 
candidate [PFH2005] 295, which is discussed later 

We found 149 sources with a significance for variability > 
3.0 out of the 300 examined sources. There is a bias towards 
bright variable sources, because for bright sources it is much 
easier to detect variability than for faint sources. 

Table |6] lists all sources with a variability factor larger than 
five in descending order The source number (Col. 1), source 
name (2), maxima of flux variability (3) and maxima of the 
significance parameter (4) are given corresponding to Table 5 
(Cols. 1, 152, 148 and 147). The next column (5) indicates the 
type of the source. If Fvar > 10, sources formerly classified as 
<hard> are now classified as <XRB>. Time variability can also 
be helpful to distinguish between foreground star and SNR can- 
didates. In some cases we had to change the source type with re- 
spect to PFH2005. This is indicated in the comment column (6). 
Column 6 also contains references to the individual sources in 
the literature. In some cases the reference provides information 
on the temporal behaviour and a more precise type (see brack- 
ets). The numbers given in connection with IVoss & Gilfanovl 
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Figure 4. Variability factor of M 31 centre sources of PFH2005 and Sect. 3 in the 0.2^.5 keV band comparing average fluxes of 
the XMM-Newton EPIC observations from June 2000 to July 2004 plotted versus HRl in the left panel and HR2 in the right panel, 
respectively. For source classification see Fig.|3] Sources with a statistical significance of the variability below 3 are marked in green 
(grey). 
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Figure 3. Variability factor of M 3 1 centre sources of PFH2005 
and Sect. 3 in the 0.2^.5 keV band derived from average fluxes 
of the XMM-Newton EPIC observations from June 2000 to 
July 2004 plotted versus maximum detected flux (erg cm"^ s"'). 
Source classification from PFH2005 is indicated: Foreground 
stars and candidates are marked as big and small stars, AGN can- 
didates as small crosses, SSS candidates as triangles, SNR and 
candidates as big and small hexagons, GlCs and XRBs and can- 
didates as big and small squares. Sources with a statistical sig- 
nificance for the variability below 3 are marked in green (grey). 
Source numbers from PFH2005 and Sect. 3 are indicated for 
sources with flux variability above 5 or maximum flux above 
8x10"'-^ erg cm"^ s"'. 



(I2007h and lWiUiams etaH (l2006bh are the Chandra derived vari- 
ability factors obtained in these papers. From the 44 sources 
listed in Table|6] six show a flux variability larger than 100. With 
a flux variability factor > 830 source 883 is the most variable 
source in our sample. Source 335 has the largest significance of 



variability, with a value of « 85. Only for ten sources the signif- 
icance of variability is below 10, for two below 5. Twenty-eight 
sources are XRBs or XRB candidates and seven are SSS candi- 
dates. 

Table |7] lists all "bright" sources with maximum flux larger 
than 8x10"'"^ erg cm"^ s"' and a flux variability smaller than 
five, giving the same information as in Table|6] All seven sources 
listed in Table Q have a significance of variability > 10. Apart 
from source 341, they are XRBs (three in globular clusters) or 
XRB candidates. The most luminous source in our sample is 
source 297 with a luminosity of ^ 3.6x10^*^ erg s"'. 

Figure |4] shows the relationship between the variability fac- 
tor and the hardness ratios HRl and HR2, respectively. We used 
the hardness ratios of the observation from which the source en- 
tered the catalogue of variable sources. The HRl plot shows that 
the sample of highly variable sources includes SSS and XRB 
candidates, which occupy two distinct regions in this plot (see 
also Haberl & Pietsch 1999, for the LMC). The SSSs marked 
by triangles, appear on the left hand side, while the XRBs or 
XRB candidates have much harder spectra, in agreement with 
their classification. In the HR2 plot the highly variable XRBs 
and XRB candidates are, apart from the two sources classified 
as <SNR>, separated from the bulk of the less variable sources 
by sources classified as <hard>. Due to the distinct temporal 
variability of these sources and the strong absorption in the cen- 
tral region of M 31, it is very unlikely that they are AGNs. So 
only <fg star> or <XRB> will be left as possible classification. 
In accordance with the hardness ratios we suggest sources 169, 
225, 322, 328, 335 and 420 as XRB candidates. 

Individual sources are discussed in the next section. 



5. Discussion 

In each of the following subsections, we first discuss the sources 
described in the catalogue extension (Sect. |3]l. In addition we 
reclassified some sources of PFH2005 based on the results of 
our time variability study and on recent papers in the literature. 

We classified the sources described in the catalogue exten- 
sion into different types of X-ray emitting objects: foreground 
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Table 6. Variable sources with flux variability larger than 5, ordered by variability. 
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Source 


Name 
XMMM31 J 


fvar 


svar 


fmax'' 


type+ 


Coniment^ 




004947 S-i-Al 111^ 

UUH-Z.H- / . O T^H- 1 i 1 i 


811 10 

O.^ 1 . IVJ 


S4 7S 


18 09 




1 {v\ 1(\ 09 9 19 17 


jyyj 


OnA^O'> 7j-A1 1 70^ 


uZ'+.Uj 


70 81 


49 71 
4Z. / 1 


<- A. rvD -> 


1 ^'t^ 9/"! 054 9^ l*"! 91A1^ 90 l^^t^ 91 

z(^i, yj4.zj, ziojj, zu, zii^ij, zj 


787 


0049^4 ^+41 1 SIO 


isi fn 

J J J.U / 


41 9A 


91 OS 


/v fv D 


ll't'i 90 170 S'^ ISO RT4 VRN'i 910^ 99r'vt1 


J i\j 


004949 1 4-41 1 ^^08 


901 71 

Z.U 1 . / 1 


A9 44 


88 47 


yVl\.D 


](t'\ 90 4A8 S'^ lO 985^ lS('\7t RH VRT^ 1 OO^i 99('vt'i 


o / o 


004 14A 7a.41 1111 
UU'T-lH-'+. / T'+l 1111 


1 78 7Q 
L / o. /y 


40 90 


IS Al 
J J. 01 




l(,t,sv;, 4(,t; 




004^00 8-4-41 1 QOO 


1 11 71 


^^7 99 
J / .zz 


14 9 S 


^VRR-i 
^zvlxO 


^AmV--*^ 1(\ OA 1^ li't 107^ 10 M(\7\ 11 \A(\7\ 90 
1 <-/\Oi> > }, Z(^L, yo. J }, J 1,1, LxJ / ), lU, 1 1 J, 14^^V Z,\J, Z,Z,\\,lj 


909 


00490S 84-41 1 ^90 


07 80 
y 1 .oy 


99 41 


1 1 19 
U. JZ 


^YDR^ 
<^yVl\.D^ 


](r\ 1(\ 1C\9.\ lO 01^ 19 1 Si't^ 910^ 




004^^07 1 4-4 11810 

UUH-^U / . 1 tH- 1 1 O 1 U 


07 fx'> 
y 1 .OJ 


9S 97 
ZJ.Z / 


1 7S 
lU. IJ 


^VRR-^ 
<^A.1S.1j> 


1(\ Afx \ \ l('t 1 SS^ 90 94 
H^tJ, Z(,l, 40. 1 J, Jl,l, 1 J Jj, ZU, ZH^l^, Z4 




004T18 8-1-419017 
UU'4-^lo.o+'4-lZUl / 


8^ Ql 


40 S4 


1 A IS 
14. JJ 


^CCC^ 


1 (r\ 9 l('t OA^ 1 11 1 A(^7\ 1 S^'lA 1C{{\t\ 99('\A 

\-\X)^ z, j(^i, yoj, lu, ij, i4(^v^, ijt^v^, zui^v^, zzt^v^ 


A.'X 1 

H- J 1 


004^1 S-i-41 1 l'^f\ 


89 f\9. 


40 17 

H-U. 1 / 


1 7A 
1 y . / u 




lO^ lO A04'^ ISi'vt'i 


ool 


004941 84-41 l^^^'^ 


7A 97 
/ u.z / 


70 A8 

/y.oo 


8A 1 7 
oO. 1 / 


VRR 


H^L,i,SV ^, 4^^L^, O^^L, LiVlAD }, lU, ZZ^^V J 


OJD 


0049 74-41 1 84^^ 


^^7 41 


1 9 AO 

iz.oy 


8 19 
o. JZ 


^VRR-^ 
<^A.lS.D-> 


l(t 70 0^ l/'t 9A0^ 1 S/'t^ 10 91('t^ 99('\7t'\ 

zi^L, /y.uj, j^i, zou^, iJi,i^, ly, zn^ij, zzi^v,i^ 


1 AO 


004 1 4^ 19118 


44 A7 


90 99 
ZU.ZZ 


77 
y. / / 


<AK_D> 


former type; <hard>; 1, 14, 24 


SS7 


0049 S9 4j-41 1 ^^4Q 
UU''tZJZ.'+t''+1 LKy-ty 


IQ AO 


90 98 
ZU.Zo 


10 

y.oy 


^VRR-i 
^yvlxD> 


Ift A4 A^ 
Z.\i, 04. O^ 


■^90 


00494'^ 1-1.41 179T 
UU'4-Z'+ J . 1 -TH- 11/ Z^ 


IS 07 


94 94 
Z'+.Z^ 


1 1 71 
1 1 . / 1 


<Alx£S> 


1 fr c^7\ 9/'t OQ S^ l/'t 1 SS'l 99('\/ 


J 1 o 


00494^ ^4.41 1 ^1 Q 
UUH-ZH- J . J T'+ 1 1 J 1 7 


14 47 


9 1 10 
Zl . JU 


A 17 
0.1/ 




loiiiier lype. <.oiMxP', z^^l^, zu, zz, Z4 


J J J 


004947 1-1-41 1A9Q 


14 09 


84 AO 

04. oy 


1 4A 44 
140.44 


<A.Jv£S> 


fnrrrK^r t\Trt(^- ^Vinr-r\~~^- 1 /"ciA 9 10 19 11 14 90 99('\A 

lormer type. <narQ>, ii^sv^, z, lu, iz, u, 14, zu, zzt^v^ 


ooo 


004^00 Qj-41 9^^9 


10 SO 


1 8 9A 
lo.ZO 


A 7S 
4. / J 


VRR 


1 (t\ 9/'t'\ 4/^1^ 7('t T A/TVR^ 
H^L^, Z(,L^, 4^^L_f, /^^L, i-,l\iJ\.D} 


o / J 


004T 1 8 7 J.4 1 1 804 


98 78 
Zo. / o 


1 O 70 


A S7 
4. J / 






ooU 


0049^^ 04-419^^1 


98 09 
Zo.UZ 


1 10 

lu. jy 


9 A8 
Z.Oo 


^VRR-^ 
<^A.1\.D> 


l(t AS 9^ 
Zl^l, OJ.Z } 




004T 1 4-l4 1 1AAC\ 


97 OA 
Z / .UD 


90 07 
ZU.U / 


S 14 
J. J4 


YRR 




S7Q 

o ly 


004994 '^j-4 19401 


94 Q'\ 


Q 1^ 

y.lj 


1 09 

1 .yz 


^YRR-i 
^zv1x1j> 


z 




0049^^^ 1j-A \ 1 ^40 


9 1 Q7 

LV.y 1 


1 1 77 
11.// 


S SA 

J. jO 


<igijiai> 


f/^rm*ir t\7t-**i- ^QMR->.- 9 1114 90 99 94 

loimei type. <oi>K>, z, u, 14, zu, zz, Z4 


99^^ 


0049 1 Q-i-4 1 C\f\Al 


1 S 7A 

1 J. /D 


8 10 


S 90 

J.ZXj 


^YRR-^ 
<^A.1S.1j> 


former type; <hard>i 2, 22(v) 


S7 1 
o / 1 


0049^4 A-x,4 11819 
UU'4-Zj4.0t'41 lolZ 


1 S OO 


1 A7 
J.O / 


1 ss 

1 .J J 


^ccc^ 


1 8 
lo 


1 Ql 

VyV 


0041 S4 ^-1-410794 


14 18 
1'+. 1 o 


14 01 


A 48 
4.4o 




1 (t\ 


AOC\ 
H-ZU 


004^ 1 fx 0-1-4 1 1 849 
yjXJ'-tJ 1 D.Ut'+I 1 OH-Z 


1101 
1 j.Ul 


1 7 

1 / .J J 


7 00 

/ .yy 


^YRR-i 
<- A. rvD ^ 


luriiiei type, ■^iidiu^, i, z, ij, zu, z,z,\\) 


^99 


004944 94.41 9 80Q 
UU'+-Z'+'+ .Zt'+IZoU" 


1 9 OA 


1 1 81 


A 14 
4. J4 


<-y\ fv D 


former type; <hard>i 1, 2, 13, 14 


^98 


00494^^ 04-41 1407 


I 1 AO 

I I .oy 


1 Q 99 
ly.zz 


A 90 

o.zy 


^YRR-^ 
<^A.1S.D> 


fr\rmfir f^^n^i- ^^ViQrH^- 1 9 19 11 90 99 

loriiier type. <.iiai(j>, i, z, iz ,ij, zu, zz 


^90 


00494^ 84-41 17SA 
UUH-Z'+j.ot'+I 1 / JO 


1 70 
lU. /U 


1 1 80 

1 j.oy 


1 S8 

1 .JO 




l*"! S1^ 90 94 
Jyi, J i^, ZU , Z4 


oO'+ 


004947 Q-l41 1 ^^40 

UU'4-Z'+ / .y+'+i 1 j'+y 


1 O 94 


9S 8'> 


41 

y.4i 




9 90 99 94 

Z, ZU, ZZ, Z4 




004^^08 ^^-1-41 1 890 
UUH-JUo. JT'+l 1 oZU 


1 O'^ 


1 9 08 
IZ.Uo 


1 09 

1 .yz 


^^^^^^ 


If'f 18^ 
J\l, Jo J 


SS9 
OOZ 


004949 0-1-41 1 '^.'X'X 
UU'4-Z'=+Z.U+'4-l 1 J 


Q 9A 

y.zo 


1 4 18 
14. JO 


1 o so 

lU. JU 




9 90 99('\A 
Z, ZU, ZZ(^V ) 




004^^^^ ^^j-41 081 9 


8 97 
o.z / 


8 S4 
o. J4 


9 A8 
Z.Oo 


oil- 


J 


9^^ 


004991 ^^-1-41 1418 
UU'^-ZZl.O+'+l I'+lo 


7 48 


1 7 SA 
1 / .00 


7 11 
/ . J 1 




1 Knrct^ 9 8 90 99/^l7^ 

H^sv,Dursi^, z, o, ZU , zz{y) 


9^0 


004919 1-1-41 17S7 

UUt^Z. 1Z,.1Tt'1 1 1 J 1 


A 45 


27.71 


21.25 


<GIC> 


Wwl 9 S 19 1 S('v1 lA 90 

Lyciv Z,, J, IZ, 1J^V^,1U, ZU, ZZ^^ V 


316 


004242.8+411639 


6.37 


8.01 


4.00 




former type; <SNR>; 10, 12(v), 20, 22(v) 


384 


004303.3+411527 


6.19 


37.24 


35.54 


<XRB> 


l(sv), 2(t, 58.6), 3(t, 33), 5, 10, 12(v), 13, 14, 20, 22(v,t) 


465 


004333.4+412140 


5.58 


5.13 


2.28 


<hard> 


2 


865 


004323.4+412208 


5.52 


3.57 


0.63 


<SNR> 




208 


004207.0+410017 


5.35 


5.10 


8.75 


<GIC> 


5, 13, 14, 16, 21 


249 


004219.6+412153 


5.35 


10.47 


1.51 


GIC 


2, 5, 16, 22 


415 


004314.5+411649 


5.33 


5.21 


1.03 




2, 22, 24 


858 


004250.4+411556 


5.14 


9.69 


3.10 


SNR 


2, 9, 22, 24 



Notes: 

* : maximum flux in units of Ix lO"''' erg cm"- s"' or maximum luminosity in units of 7.3x 10^^' erg s"' 

* : type according to Table[2l partly changed as mentioned in the comment column 

t : 1: U rudolyubov et al.' ('20061^, 2: 'Voss & Gilfanov llOOf), 3: 'Williams et al.' /2006B), 4: 'Trudolyubov et al.' ('200 6j), 5 : 
iTrudolvubo v & Priedhorskv (2004), 6: Williams et al. (2006a), 7: Williams et al. (2005b), 8: Pietsch & Haberl (2005), 9: Kong et al. ( 2003bh . 
10: iTrinchieri & Fabbianoi (l99lh, 11: ICollura et alj ^1990), 12: IPrimini et al.l a993l), 1 3: ISupperetalJ (119971) 14: ISupper et al. (200T]), 15: 
lOsbome et aD l l200ll) , 16: lDi Stefano et all ll200i), 17:lFan et alj j2005l), 18: lPietsch et aP ll2007h, 19:lGarcia et aP j200(]l) , 20: fKaaret l2003) . 21: 
IWilliams et all ( l2004h . 22: iKong et al.l ( |2002|) . 23: IWilliams et all j2005al) . 24: iDi Stefano etal] ( |2004|) . 25: iBamard et all l l2003al) : t: transient, v: 
variable, sv: spectrally variable, r: recurrent, d: dipping, z: Z-source candidate; BH: black hole, XRN: X-ray nova, XRT: X-ray transient, LMXB: 
low mass X-ray binary, NS: neutron star; numbers indicate the variability given by the corresponding paper 



stars (fg Star), galaxies (Gal), AGN, supersoft sources (SSS), 
supernova remnants (SNR) and X-ray binaries (XRB), using the 
X-ray properties together with information from catalogues at 
other wavelengths. The selection criteria for these classes are 
given in Table ID Additionally we use the time variability to 



classify sources. In the field of M 31 mainly XRBs or SSSs 
can show very strong variability (Fvar ^ 10) on time scales 
of years. In only a few cases we were able to identify an X- 
ray source with a source already classified from the optical, in- 
frared or radio data. We have no well-defined hardness ratio cri- 
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Table 7. Sources with maximum flux larger than 8x10 erg cm ^ s a statistical significance of variability larger than 10 and a 
flux variability smaller than 5, ordered by flux. 



Source 


Name 


fvar 


svar 


fmax* 


type^ 


Comment^ 




XMMM31 J 












297 


004238.5+411603 


1.56 


47.20 


49.71 


XRB 


l(sv,z), 2, 10(v), 12(v), 13, 14, 20, 22(v), 25(LMXB) 


257 


004223.0+411534 


3.05 


51.35 


16.84 


<XRB> 


l(sv), 2, 10(v), 12(v), 13, 14, 20(v), 22(v) 


239 


004215.7+410115 


1.48 


10.70 


16.73 


GIC 


10, ll(v), 12, 13, 14(v), 16, 20, 21(v) 


408 


004310.6+411451 


1.37 


12.81 


10.77 


GIC 


l(sv), 2, 5, 10, 12, 13, 14, 16, 20, 22(v) 


353 


004252.5+411854 


2.07 


29.97 


9.68 


<XRB> 


l(sv), 2, 10, 12, 13, 14, 20(v, NS-LMXB), 22(v) 


341 


004248.5+411522 


1.27 


10.99 


8.94 


<hard> 


l(sv), 2, 10, 12, 13, 14, 20, 22(v,sv) 


414 


004314.3+410722 


2.47 


26.89 


8.21 


GIC 


l(d,sv), 2(t, 53.4), 5, 10, 12, 13, 14, 16, 20, 22 



Notes: 

* : maximum flux in units of 1 x 10"''' erg cm"^ s"' or maximum luminosity in units of 7.3x 10^* erg s"' 

* : type according to Table[2] 

' : for comment column see Tabled 



teria to differentiate between <hard> sources (XRBs, Crab-like 
SNRs or AGNs). Fifteen sources of the catalogue extension are 
classified as < hard> (see Table [ 2 ]l. Three of them were found 
with Chandra jKong et alj|2002l: IVoss & Gilfanovll2007h . Five 
sources remain unidentified or without classification. Two of the 
five are already kn own from Chandra observations (see Table|6]). 
iKongetalJ ( 120021) classified source 884 as SSS. 



5.1. Foreground Stars 

Foreground stars are a class of X-ray sources which is homoge- 
neously distributed over the field of M 3 1 . The good positioning 
of XMM-Newton and the available catalogues USNO-B 1 and 
2MASS allow us to effectively select this type of source. We 
found one foreground star candidate (877) in our source cata- 
logue extension. From the optical colours in the USNO-B 1 cata- 
logue we estimate the typ e to be A3 III o r A5 III, using the stellar 
spectral flux library from lPicklesI (Il998l) . Another possible fore- 
ground star candidate (859) is a USNO-B 1 and 2MASS source. 
From the USNO-B 1 magnitudes we derived /x //opt ~ -0.87 and 
/x//opt,R ~ -1.27, where /opt.R is the flux in the R-band. The 
/x//opt value is too large, to satisfy our foreground star selec- 
tion criterion. But for very red objects it can be sufficient that 
/x//opt,R is < -1. The source could be a foreground star, in 
agr eement with the va lues we found for source 295 (see below). 
But iKim et akl (|2007|) suggested this optical source as possible 
globular cluster. This classification would also be in agreement 
with our hardness ratios, fx/fopt values and USNO-Bl magni- 
tudes (see [pan et al. 2005i). So we cannot decide on a fg Star or 
XRB nature and we classify source 859 as <hard>. 

PFH2005 classified source 295 as a SNR. This classification 
has to be rejected due to the distinct time variability we found. 
We created light curves in the 0.2 - 2.0 keV range for the differ- 
ent observations. In some, especially in c3 (see Fig.|5]l and in c4, 
the source shows strong flares. The observation c2 consists of the 
decaying wing of a strong flare, while the source remains rather 
quiet in "b". In addition we carefully chec ked the 2MASS and 
Local Group (LG) survey R-band images (iMassev et al.ll2006l) 
and found in both images a faint point-like source, at the X-ray 
position. Eq. © gives /x//opt ~ -0.66 and /x//opt,R ~ -1-28 
using brightnesses from the LG survey photometric catalogue. 
The /x //oDt values derived from the catalogue by Haiman et alj 
(119941 /x//opt ~ -0.60 and /x//opt,R ~ - 1 .44) are in good agree- 
ment with the values derived form the LG Survey and are reason- 



Time (s) 



Figure 5. Summed EPIC PN, MOS 1 and MOS 2 0.2 - 2.0 keV 
light curve of source 295 in the c3 observation binned with 1000s 
and without background subtraction. The zero time corresponds 
to 2001-06-29 07:53:36. 



able for a red star. Considering all those points, we now classify 
this source as a foreground star. 



5.2. Supersoft sources 

Spectra of SSSs with low energy resolution can be modelled by 
black body spectra with temperatures below 50 eV. They radi- 
ate close to the Eddington luminosity of a 1 Mq object and are 
believed to be white dwarf systems steadily burning hydrogen at 
the surface. They were identified as a class of X-ray sources by 
ROSAT and a re often observed as transient X-ray sources (see 
lGreineijl2000i and references therein). 

Our catalogue extension contains three SSSs. Two of them 
(871, 886) correlate with optical novae and have been investi- 
gated in more detail in PFF2005 and PHS2007. 

The third one (875) lies near source [PFH2005] 431 (dis- 
tance X 12"). As source 431 is brightest in observation cl and 
source 875 is detected in observation c4, we can exclude that 
they are the same source. From the time variability and the posi- 
tional errors it would be possible that source 875 corresponds to 
the nova M31N1923-12b (= [H29] N28; distance « 7"), which 
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was reported in the optical wavelength regime bv lHubbld([T929l 
see also Nova list of PHS2007). However super soft X-ray emis- 
sion from novae up to now has only been observed up to ten 
years after the optical outburst (see e. g. PHS2007). So if source 
875 really coincides with M31N1923-12b, the X-ray emission 
we found would have to be connected with an unreported opti- 
cal outburst, which occurred during the last ten years, making 
the source a recurrent nova. Another possibility is, that source 
875 corresponds to yet another nova, not reported in the optical. 
But we cannot exclude that source 875 is not a nova at all. 

5.3. Supernova remnants 

SNRs can be separated into sources where thermal components 
dominate the X-ray spectrum below 2 keV, and so-called "pleri- 
ons" with power law spectra. The former are located in areas of 
the X-ray colour/colour diagrams which only overlap with fore- 
ground stars. If we assume that we have identified all foreground 
star candidates from the optical correlation and inspection of the 
optical images, the remaining sources can be classified as SNR 
candidates using the criteria given in Table |2] 

We thus identified six SNR candidates in our catalogue ex- 
tension. One of them (885 ) had been pr eviously observed with 
Chandra (iKong et alJl2002l:lKaarejl2002h . but had not been clas- 
sified. A second source (858) coincides with a source reported as 
a ring-like extended object from Chandra observations, which 
was also detected in t he optical and radi o wavelength regimes 
and identified as SNR dKong et al.ll2003bl) . 

Two sources from our catalogue extension, which are classi- 
fied as SNRs are listed in Table|6l Source 858 lies next to source 
875, which was first detected in observation b. Therefore the flux 
of source 858 is underestimated in "b" and the source appears 
variable. There is thus no need to change the type of this source. 
For source 865 we can only gain upper limits for the flux, apart 
from observation c3 (Lx ~ 4.6x10"^^ erg s"'), which leads to a 
significance of variability of only 3.57, not much above the 3 cr 
limit. So the source can still be classified as a SNR candidate, 
despite the alleged time variability. 

We now discuss the SNR candidates of PFH2005, which 
show time variability: 

Source 318 shows significant variability. Therefore we have 
to reject the classification of PFH2005 as <SNR>. Fig.|6]shows, 
that in observation "b" the source is about a factor of 10 to 35 
less luminous than in the other observations. We checked care- 
fully whether the source lies at the rim of a CCD or on a CCD 
gap. Neither is the case. In the following we discuss possible 
source classifications: the hardness ratios are in agreement with 
our foreground star criterion, however, the duration of the out- 
burst of about two years seems much too long for a stellar flare 
(Fig.|6]). Since we also di d not find an optical counterpart in the 
images of the LG survey jMassev et al.ll2006l) . we can exclude a 
foreground star identification. The behaviour on long-term time 
scales suggests an X-ray nova as a possible source c l assification 
(iHaberl & Pietschlll999l: iTanaka & Shibazakilll996t IChen et alJ 
Il997h . We used the data of observation c2, in which the source is 
most luminous, to produce an EPIC PN spectrum. A disk black- 
body model fitted to the spectrum gives a temperature at the in- 
ner edge of the accretion disk of =s 190 eV, which seems too 
small for an X-ray Nova or LMXB. We also fitted a blackbody 
spectrum. The temperature of ^ 160 eV is too high for a SSS, bu t 
would be in agreement with a QSS ( IOridl2006t lFabbiandl2006[) . 
A power law fit gives a photon index of ^ 4.7. Photon indices of 
XRBs and AGNs are much smaller than that value. So the nature 
of this source remains unclear. 



5.2x10^ 



5.25x10^ 
Modified Julian Day 



5.3x10^ 



Figure 6. EPIC long-term light curve of source 318. We used 
XID fluxes. The arrow marks a 3cr upper limit. 



For source 316 the variability factor we found is based on 
observation b (without b: Fyar = 3.07 and S'vai- = 4.75). As the 
source lies next to the bright transient source 881, which was first 
detected in that observation, the flux of source 316 may be un- 
derestimated and the source could appear as a variable. However 
due to the variability reported in the literature (see Table |6]l, the 
SNR classification has to be rejected. 



5.4. Globular cluster sources and X-ray binaries 

A significant part of the luminous X-ray sources in the Galaxy 
and M 31 are found in globular clusters. We correlated our cata- 
logue extension with that of Galleti et al. (2004). 

All <hard> sources of our source catalogue extension, which 
have a variability factor larger than ten are classified as XRBs. 
References for these sources can be found in Table |6] TPC06 
report on four bright X-ray transients, which they detected in 
the observations of July 2004 and suggest them as XRB candi- 
dates. We also found these sources and classified source 878 and 
identified sources 881, 888, 890 as XR Bs. One of th e identified 
XRBs (890) shows a very soft spectrum. IWilUams et al. (2005^ 
observed source 888 with Chandra and HST. From the location 
and X-ray spectrum they suggest it is a LMXB. They propose 
as optical counterpart a star within the X-ray error box, which 
shows a change in optical brightness (AB) of ^ I mag. Source 
881 was first detected in January 1979 by TF91 with the Einstein 
observatory. WGM06 rediscovered it in Chandra observations 
from 2004. Their coordinated HST ACS imaging does not re- 
veal any variable optical counterpart. From the X-ray spectrum 
and the lack of a bright star WGM06 suggest this source as a 
LMXB with a black hole. 

In PFH2005, sources 169, 225, 322, 328, 335 and 420 were 
classified as <hard>. We found that they all have a time variabil- 
ity factor larger than ten and therefore re-classified them as XRB 
candidates. 

Sources 257 and 384 were proposed a s stell ar mass black 
hole c andidates by iBarnard et al.l (l2003bl) and iBarnard et al.l 
( |2004|) . respectively. Recently, it was shown that the aperiodic 
variabi Uty of these sources has an artificial origin dBamard et al.l 
I2007bh . So there is no longer clear evidenc e for a black hole 
nature of these objects (IBarnard et alj|20()7ab . We now classify 
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Table 8. Outbursts of source 883 = [PFJ93] 51 = 


rVG20071 136 


— [ Ir^^/UUOJ / / 








Satellite 


Time of observation 




Reference' 


ROSAT HRl 


Jul 1990 


4.7 


1 


ROSAT HRl 


Jul/ Aug 1995 


4.6 


2 


Chandra AClS-1 


May 2004 


3.3 


3 


XMM-Newton EPIC 


Jul 2004 


3.7 


4, this work 



Notes: 

* : 0.5 - 8.0 keV unabsorbed luminosity in units of 10^' erg s"' for 
a distance of 780 kpc, assuming A'h = 1.1 xlO"' cm"^ and a photon 
index of 1.6 

t : 1: iPrimini et a □ (Il993h, 2: IRX H catalogue, 3: IVoss & Gilfaiiovl 
( 120071) , 4: lTnidolvubov et alj ( l2006bl) 



sources 257 and 384 as XRB candidates, based on their time 
variability (see Tables |6] and |7]i. 

Source 883 is a transient, only detected in July 2004 (obs. b) 
in our study. It stands out in Fig.[3]and Table|6]as the source with 
the highest variability (Fyai ~ 830). The EPIC pn data of source 
883 during observation b can be well fitted with an absorbed 
power law model (A^h = 1.1 + 0.2 xlO~' cm"-, photon power 
law index = 1 .61 + 0.08 , unabsorbed 0.5 - 8.0 keV luminosity - 
3.7x10^' erg s ~'). The source cor relates with the GIC candidate 
Bo 128 (e.g. iGalleti et ai]|2004 . Based o n its variabiUty, lu- 
minosity and absorbed power law spectrum iTrudolvubov et all 
(2006b) classify the source as a neutron star XRB candidate 
(#77 in their list of bright X-ray sources detected in the cen- 
tral part of M 31). During the Chandra monitoring of the centre 
area of M 3 1 the transient was detected at a si milar luminosity 2 
month s earlier in May 2004 (source 136 in IVoss & Gilfanovl 
l2007b . most likely during the same outburst. No additional 
Chandra detections of the source have been reported. No source 
was detected at the position of this bright transient with the 
Einstein Observatory 1979/80 (e.g. TF91), during the ROSAT 
PSPC sur veys (Jul 2001, Jul / Aug 2 002, Dec 2002/Jan 2003, Jul 
2003; see lSupper et'aLlll997ll200lh and during ROSAT HRl ob- 
servations in Jul 2004 and Jan 2006 (see source catalogues of 
the pointed HRl observations IRXH). However, two additional 
outbursts of the transient were detected with the ROSAT HRl in 
July 1990 (source 51 inPFJ93) and in Jul/Aug 1995 (see IRXH). 
The luminosity derived for these outbursts is remarkably similar 
to the luminosity of the outburst in 2004 if we assume that the X- 
ray spectrum of this recurrent transient can always be described 
by the same model as during the 2004 outburst (see Table[8]l. 

6. Conclusion 

In this paper we present an updated source list of the central area 
of the bright Local Group spiral galaxy M 31, using the obser- 
vations from June 2000 to July 2004 available from the XMM- 
Newton archive. We extended the source catalogue by PFH2005, 
based on the merged images of the observations from 2000 to 
2002 by searching sources in the observations of 2004 and re- 
examining the observations used in PFH2005 individually. To 
classify or identify more of the sources, we examined their long 
term time variability. 

We obtained 39 sources in addition to the 265 reported by 
PFH2005 in the field. The identification and classification of 
these sources is based on properties in the X-ray wavelength 



regime: hardness ratios and temporal variability. In addition, in- 
formation from cross correlations with M 3 1 catalogues in the 
radio, infra-red, optical and X-ray wavelength regimes are used. 

We detected three SSS candidates, one SNR and six SNR 
candidates, one GIC candidate, three XRBs and four XRB can- 
didates. Additionally we identified one foreground star candi- 
date and classified fifteen sources as <hard>, which may either 
be XRBs or Crab-like SNRs in M 3 1 or background AGNs. The 
remaining five sources remain unidentified and without classi- 
fication. Two sources were found to be extended. One of them 
was classified as <hard>. The other stays without classification. 

To examine the time variability we calculated the flux or at 
least an upper limit at the source positions in each observation. 
We determined the variability factor and significance parameter 
for each source, comparing the XID flux ratios of the different 
observations with each other. The time variability helped us to 
decide if a source classified as <hard> in PFH2005 can be an 
XRB candidate. In addition we could use time variability to dis- 
tinguish between foreground star and SNR candidates. 

Six sources of PFH2005, which were classified as <hard>, 
show distinct time variability. Based on that variability, their 
hardness ratios and the strong absorption in the centre of M 3 1 
we suggest these sources as XRB candidates. The SNR classi- 
fication from source 295 was changed to foreground star due to 
the distinct time variability we found and its identification with a 
faint stellar object. Other SNR classifications (sources 316, 318) 
were rejected due to time variability of the sources. 

To verify our suggested classifications further investigations, 
including at other wavelengths will be necessary. 
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